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ABSTRACT 

We present numerical testing results of our mass function derived in our 
previous paper, and compare the testing results with those of the popular 
Press-Schechter (PS) mass function. Two fiducial models are considered for 
the test: the scale-free power-law spectra P(k) oc k n with spectral indices 
n = —1,0 and the standard cold dark matter (SCDM) model with Q = 1 and 
h = 0.5. For the power-law models, we use numerical data averaged over several 
different output times: ten output times from two N-body realizations for the 
n — — 1 power-law model; four outputs from one realization for the n = model. 
While for the SCDM model, we consider four outputs separately at redshifts 
z = 0,0.43, 1.14 and 1.86 from one large N-body simulations. The comparison 
results show that our mass function fits the numerical data in a much improved 
way over the PS one. Thus, we expect that our mass function can be a viable 
alternative of the PS mass function in applications to various areas. 

Subject headings: cosmology: theory — large-scale structure of 
universe 
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1. INTRODUCTION 

One of the most fundamental problems in cosmology is to understand the formation 
and evolution of the large-scale structure in the universe such as galaxies, groups and 
clusters of galaxies, etc. In order to understand the large-scale structure, however, it 
is highly desirable to have an analytical framework within which theoretical predictions 
for structure formation can be made. The cosmological mass function, n(M) [n(M)dM: 
defined as the comoving number density of gravitationally bound structures - dark halos 
with mass M\ provides this analytical tool since different candidate models for structure 
formation predict different number densities of dark halos. 

Press & Schechter (1974, hereafter PS) developed for the first time an analytic 
formalism to evaluate the mass function. Finding a mass function requires both dynamics 
and statistics. Dynamically PS adopted the top-hat spherical model, according to which the 
collapse condition for forming dark halos is determined purely by its local average density. 
Statistically PS assumed that the initial density field is Gaussian, and selected dark halos 
from the peaks of the linear Gaussian density field. 



The practical success of the PS mass function (e.g., [Efstathiou et al. 1988| ; [Lacey fc 



Cole 1994| ) and the absence of viable alternatives led many authors to use it routinely in 



the last decade(e.g., Cole fc Kaiser 19881 ; [White fc Frenk 1991fc Kauffmann et al. 1993| ; [Fan 



Bahcall, fc Gen 1997|) . 



But, recently high resolution N-body simulations have shown the limitation of the 
PS theory: First, it has been clearly shown by several N-body simulations that the true 
gravitational collapse process must be nonspherical (e.g., [Shandarin et al. 1995| ; [Kuhlmanl 



et al. 1996|) , indicating the weakness of the PS dynamical background. Second, N-body 



simulations also showed that peaks of the linear density field are "poorly" correlated with 
the final spatial locations of dark halos (e.g., [Katz, Quinn, fc Gelb 1993]) , although the PS 
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formalism assumes that dark halos form in the peaks of the density field. Third, recent 
high-resolution numerical tests have detected that experimental results are flatter than 
the standard PS mass function in shape (e.g., Governato et al. 1998; Tormen 1998 and 
references therein ; Bheth fc Tormen 19991 ) . 



Currently various attempts have been made to find a better-fitting mass function. 
One approach to better mass functions has been focused on finding phenomenological 
fitting parameters, keeping the original PS formula unchanged; for example, regarding the 
density threshold as a function of redshift, etc (see Governato et al. 1998; |Sheth fc Tormen 



|1999[ ). Another approach has been on improving the dynamics of the PS formalism by 
implementing anisotropic collapse conditions (e.g., Monaco 1995| ; Audit, Teyssier, & Alimi 



1997). A full analytical alternative of the PS mass function using this approach has been 
found by |Lee fc Shandarin ( 1998 )| , in the frame of the Zel'dovich approximation. This 



approach goes along with the assumption that dark halos correspond to the third axis 
collapse. Indeed the mass function we found in the previous paper is both a dynamically 
and statistically improved version of the PS one. First, it is based on a more realistic 
nonspherical dynamical model. Second, the underlying statistical assumption that dark 
halos form in the local maxima of the smallest eigenvalue of the deformation tensor, 
A3 (see §2) is in general agreement with the N-body results performed by |Shandarin fc 



Klypin (1984)| . Third, and most importantly our mass function was shown to have desired 



properties like a lower peak and more high-mass halos. That is, our mass function is flatter 
than the PS one. 

In this Letter we present numerical testing results of our mass function for the case 
of two fiducial models: the scale-free power-law spectra with spectral indices n = —1,0, 
and the standard cold dark matter (SCDM) model with Q = 1 and h = 0.5. In §2 we 
briefly summarize the analytic mass function theories for the readers convenience . In §3 we 
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explain the N-body simulations used to produce the numerical mass functions, and compare 
the analytical mass functions with the numerical results. In §4 we draw a final conclusion. 



2. SUMMARY OF MASS FUNCTION THEORIES 



The PS theory assumes that dark halos of mass M form hierarchically in the regions 
where the linear Gaussian density field 5 = (p — p) / p (p: mean density) filtered on mass 
scale M reaches its threshold value 5 C for collapse: 



nps(M) 



7TM 2 



d In a 



dlnM 



a 



exp 



2a 2 



The density threshold 5 C for a flat universe is originally given by the spherical top-hat 
model: 5 C ~ 1.69 (e.g., [Peebles 1993Q . But in many numerical tests it has been detected 
that lowered S c (roughly 1.5) gives a better fit in the high-mass section (e.g., |Efstathioii 
fe Rees 1988] |Carlberg fc Couchman 19891 ; |Klypin et al. 1995| ; |Bond fc Myers 1996Q . It is 
also worth mentioning that any PS-like formalism is least reliable in the low-mass section 
( [Monaco 1995| ). This numerical detection can be understood in the following dynamical 
argument: Although the top-hat spherical model predicts that the gravitational collapse to 
"infinite" density occurs when the density reaches 5 C ~ 1.69, halos in realistic case can form 
earlier by a rapid virialization process due to the growth of small-scale inhomogeneities 
( [Shapiro, lliev, fc Raga 1999| ). 



On the other hand, according to our approach (Lee & Shandarin 1998), dark 
halos of mass M form from the Lagrangian regions where the lowest eigenvalue A3 
(A 3 < A2 < Ai, 5 = Ai + A 2 + A3) of the deformation tensor dij (defined as the second 
derivative of the perturbation potential \& such that d^ = d 2 ty/dqidqj, qi is the Lagrangian 
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coordinate) reaches its threshold A3 C for collapse on the scale M: 

1 X 



riLs(M) 
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(2) 



In the original derivation of our mass function, the threshold A3 C for collapse has been 
empirically chosen to be 0.37. A similar logic used to give a dynamical explanation to 
the lowered S c of the PS formalism applies here. Although a simple extrapolation of the 
Zel'dovich approximation to nonlinear regime predicts that the formation of dark halos 
corresponding to the third axis collapse occurs at A3 C = 1, the first and the second axis 
collapse speed up the formation of halos, which would result in a lowered A3 C (see also |Audit 
|et al. 1997|) . 

In § 3, we witness that our mass function with the original suggested value of A3 = 0.37 
does agree with the numerical data quite well. 



3. NUMERICAL vs. ANALYTICAL MASS FUNCTIONS 

3.1. Comparison for Scale- Free Model 

The N-body simulations of a flat matter-dominated universe for power-law spectra 
P(k) oc k n with spectral indices n = —1 and were run by [White (1994)| using a 
Particle-Particle-Particle-Mesh code with 100 3 particles in a 256 3 grid with periodic 
boundary conditions. 

Tormen (1998, 1999) identified dark halos from the N-body simulations using a 
standard halo finder - the friends-of-friends algorithm with a linking length 0.2 [hereafter 
FOF (0.2)]. Numerical data for the n = — 1 power-law model were obtained for 10 different 
output times coming from two N-body realizations, and then a final n — — 1 numerical mass 
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function was obtained by taking an average over the 10 output values. While for the n = 
model 4 outputs from one N-body realization were averaged to produce a final numerical 
mass function. For a detailed description of the simulations, see [Tormen, Bouchet, fc White 



(1997)| . Here we use the final average numerical mass functions for comparison data. 



For the power-law spectra, the mass variance is given by the following simple form: 

where M is the characteristic mass scale [] defined by <j(M ) = 1. It is sometimes useful to 
define a filter-depending nonlinear mass scale M* related to M by M* = M (S c )~ 6 ^ n+3S) 
for a dimensionless rescaled mass variable M/M* such that a(M*) = <5 C (see [Lacey fc Cole 
TO|) . 



Figure 1 plots the fraction of mass in halos with mass M, dF/dlnM = (M 2 / p)n(M): 
our mass function with = 0.37 (solid line) against the averaged numerical data with 
Poissonian error bars, and the PS mass function with 5 C = 1.69 and 1.5 (dashed and dotted 
lines respectively) as well. The upper panel corresponds to the n = — 1 power-law model 
while the lower panel to the n = model. As one can see, our mass function fits the 
numerical data much better than the PS ones for the n = — 1 model in the high- mass 
section (M > M*). In fact |Tormen ( 1998 )| also used the spherical overdensity algorithm 



[SO (178)] as another halo finder, and showed that the numerical mass functions from FOF 
(0.2) and SO (178) are almost identical. We compared the analytical mass functions with 
his numerical data obtained from SO (178) and also found similar results. Whereas for the 
n = model, neither of our mass function and the PS one fits the numerical data well in 



x In |hee fc Shandarin (1998)|, the characteristic mass was notated by M*. But here we 



use M* to notate a slightly different mass scale. Readers should not be confused about this 
different notation of the characteristic mass. 
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the high- mass section. Yet in the low-mass section (M < M#) our mass function fits slightly 
better than the PS one for this case. 



3.2. Comparison for SCDM model 



Governato et al. (1998) provided halo catalogs produced from one large N-body 
realization (comoving box size of 500/i _1 Mpc, 47 million particles on a 360 3 grid) of SCDM 
model with Q = 1, h = 0.5 for four different epochs: z — 0, 0.43, 1.14 and 1.86 which 
are respectively normalized by cx(8/i _1 Mpc) = 1.0, 0.7, 0.467 and 0.35. They adopted the 
transfer function given by Pardeen et al. ( 1986 )| and also used the FOF (0.2) halo finder. 



For a detailed description of the simulations, see |Governato et al. ( 1998 )| . 



We derived the numerical mass functions from the catalogs by directly counting the 
number densities of halos in logarithmic scale for each epoch. In accordance with |Governato| 



et al. (1998)| , we consider halos more than 64 particles (corresponding to M > 10 14 M Q ) 
in order to avoid small-number effects of the N-body simulations. Figure 2 shows the 
comparison results. Our mass function with A3 C = 0.37 (solid line) fits the numerical data 
much better than the PS ones with S c = 1.69 and 1.5 (dashed and dotted lines respectively) 
for all chosen epochs. 



4. CONCLUSION 

We have numerically tested an analytical mass function recently derived by |Lee fc 



Shandarin (1998)1 an d compared the results with that of the standard Press-Schechter one. 



Our mass function is not just a phenomenologically obtained fitting formula but a 
new analytic formula derived through modification of the PS theory using a nonspherical 
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dynamical model. It is based on the Zel'dovich approximation taking into account the 
nonspherical nature of real gravitational collapse process, while the PS mass function is 
based on the top-hat spherical model. Consequently our mass function is characterized by 
the threshold value of the smallest eigenvalue of the deformation tensor, A3 C while the PS 
one by the density threshold, S c . 

We have shown that in the power-law model with spectral index n = — 1 and the four 
different epochs of SCDM f] our mass function with A3 C = 0.37 is significantly better than 
the PS one with 5 C = 1.69 — 1.5. It fits the numerical data well especially in the high-mass 
section (corresponding to groups and clusters of galaxies) for these two models. Furthermore 
it is worth noting that in the testing results for SCDM model our mass function agrees with 
the data well with a consistent threshold value of 0.37 at all chosen redshifts. 

On the contrary, the testing result for the n = power-law model has shown that 
there are considerable discrepancies between the analytical mass functions (both of our 
mass function and the PS one) and the numerical data in the high-mass section. The 
discrepancies with theory for the n = model, however, have been already detected QLacey 
|fc Cole 1994| ). Yet in the low-mass section our mass function fits the data slightly better for 
this case. 

Although we have tested our mass function only for two different models, given the 
promising testing results of our mass function demonstrated here, we conclude that it will 
provide a more accurate analytical tool to study structure formation. Further testings of 
the new mass function are obviously very desirable and will be reported in the following 
publications. 

2 At four different epochs of the SCDM model we effectively probe the dependence of the 
fit to the slope of the initial spectrum. 
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Fig. 1. — The square dots represent the numerical mass function with poissonian error bars. 
The solid line is the LS mass function with A3 C = 0.37 while the dashed, the dotted lines 
are the PS mass functions with 5 C = 1.69, 1.5 respectively. The upper and the lower panels 
correspond to the n = —1 and the n = power-law models respectively. See also the top 
left panel of Fig. 2 in |Tormen (1998) . 
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Fig. 2. — The square dots represent the numerical data for the case of SCDM model with 
n — l,h — 0.5. The solid line is our mass function with = 0.37, and the dashed, the 
dotted lines are the PS mass functions with 5 C = 1.69, 1.5 respectively. 



